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We theoretically investigate the heating rate and spin flip lifetimes due to near field noise for 
atoms trapped close to layered superconducting structures. In particular, we compare the case of a 
gold layer deposited above a superconductor with the case of a bare superconductor. We study a 
niobium-based and a YBa2Cu307_x (YBCO)-based chip. For both niobium and YBCO chips at a 
temperature of 4.2 K, we find that the deposition of the gold layer can have a significant impact on 
the heating rate and spin flip lifetime, as a result of the increase of the near field noise. At a chip 
temperature of 77 K, this effect is less pronounced for the YBCO chip. 
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I. Introduction 

In the area of magnetic trapping of ultracold atoms, 
considerable attention has been recently devoted to the 
interaction of atomic clouds with the surfaces of both 
superconducting atom chips [H, d d, i, H, H 0, H, 
M. [lOl. [III. ITa and superconducting solid state devices 
[131 U la, IB] • Technological advances will allow a new 
generation of fundamental experiments and applications 
involving the control of the interface between atomic 
systems and quantum solid state devices. The imple- 
mentation of such technologies depends on the ability to 
control and efficiently manipulate atoms close to super- 
conducting surfaces. However, below a certain separa- 
tion the atom-surface coupling is strong enough that the 
trapping potential is modified by the interaction of the 
atomic magnetic moment with the near-field magnetic 
noise. This leads to atom heating and spin flip induced 
atomic loss, which is the basic limitation in metal-based 
atom chips [13, [S [H, US The origin of near field 
noise lies in the thermally-induced fluctuating currents 
and in the resistivity of the chip materials in accordance 
with the fluctuation-dissipation theorem. The magnetic 
field noise is significantly smaller in the vicinity of a su- 
perconductor 0, 3 , however for the small distances in- 
volved, heating and thermally-induced spin flip may be- 
come relevant. 

The reduction of the spin flip rate has been shown ex- 
perimentally for different chip types: with and without 
a gold layer above the superconductor [l^, UIl, • The 
purpose of the top gold layer in atom chips is to use them 
in mirror-MOTs This can be achieved by coating 

the chip surface with a reflective metal film. The metal 
film can also protect the superconductor when it is op- 
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erated near the critical current density, to avoid possible 
quenching 7]. However, the presence of the metal film 
above the superconducting layer may increase the near 
field noise significantly. 



(a) Bare Superconductor chip structure 




(b) Gold + Superconductor chip structure 




FIG. 1: Schematic representation of the chip structure: (a) 
the bare superconductor chip structure consists of a layer of 
superconducting material of thickness d deposited on a chip 
base (usually a crystal); (b) the superconductor-|-gold chip 
structure consists of a thin film of gold of thickness di on 
top of the superconducting layer of thickness d2 < d\, which 
is deposited on the chip base. The superconducting layer is 
either niobium or YBCO. 

We report our investigation of the heating rate and 
the thermally-induced spin flip rate of a neutral atom 
held close to two different superconducting structures: 
one with and one without a metal film above the su- 
perconducting layer. We consider the case of two dif- 
ferent superconductors in the Meissner state: niobium 
(a conventional s-wave superconductor) and YBCO (a 
high temperature d-wave superconductor). We present 
a first-principles derivation of the heating rate, adop ting 
a model for the heating mechanism similar to [T^ |23| . 
For the derivation of the spin fiip rate the reader may 
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refer to earlier work [2l| . The evaluation of the rates is 
done within the framework of the electromagnetic field 
quantization in absorbing dielectric media [Ij, [1^ ■ De- 
spite the quantum electrodynamics formalism being orig- 
inally developed for dielectric and metallic media, it can 
be appropriately adopted to account for the electromag- 
netic field dissipation in superconducting materials in the 
Meissner state 15]. 

The paper is organized as follows. In Sec. [TTl we de- 
scribe the model we adopt for the heating rate and we 
give its general formulation. We present the quantum 
mechanical derivation of the heating rate in Appendix 
Rl In Sec. mil we report the formulas and simulations 
for the heating rate for two superconducting chip struc- 
tures: with and without a metal layer. The comparison 
is done by considering niobium and YBCO superconduc- 
tors. For completeness, we extend the comparison of the 
two materials by studying the spin flip lifetime in Sec. lIVI 
Conclusions are given in Sec. |Vl 



II. Model and Formulas 



The heating of an atom harmonically trapped can 
be understood as a transition from a motional state of 
the atomic trap to a higher motional state of the same 
trap. This is due to fluctuations of the trap center a;, 
as schematically represented in Fig. [21 We consider a 
three-dimensional harmonic oscillator with Hamiltonian 
H = + imcj^a;^ — a; • F where F is a force causing 
the trap centre to fluctuate. The force is proportional 
to the gradient of the magnetic potential experienced by 
the atom, which is a consequence of the Zeeman coupling 
of the atomic magnetic moment /x to a spatially varying 
magnetic field B(r). Fluctuations of the magnetic field 
B(r^) at the atom's position will induce fluctuations 
of the trap center and this process is described by the 
Hamiltonian 
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-X ■ V 



(m • B(rA)) 



(1) 



where x = Y^jXjXj, with Xj = (at + aj)/2 being 
the dimensionless quadrature operator and Xj the av- 
erage position of the trap center along the jth direction. 
In the following, a concise notation will be adopted by 
dropping the sum over j and assuming that j = x,y,z. 
For an atom initially in the state the average trap 
center for a transition to a state \m ^ n) is given by 
X = {'m\x\n) which is nonzero only for m — nil. 



hence Xj^± = yf^KjMLo^j^n -I- ^ ± i, where is the 

trapping frequency along the jth direction and M is the 
atomic mass. 

We obtain the rate Jj,n^n±i for the transition n 
nil along the jth trapping direction as derived in Ap- 
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FIG. 2: Schematic representation of the heating rate model. 
The heating rate is modelled as a transition from a given 
motional state of the trap to higher motional states. The 
trapping potential is considered harmonic. 
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(2) 



where Spi'fCAj^v.j) is the spectrum of the fluctuating 
force causing the shift of the trap centre. In order to 
obtain the heating rate, we consider the average energy 
of the system {E{t)) = P{n,t)fkj{n + i) where the 
sum is performed over the motional trap states |n), and 
P{n,t) is the probability that at time t the system is in 
the state In). The average heating rate T^j = {E{t))j is 
given by [23| 



j,n^n+l 7j,n^n — 1 ), (3) 



and by substituting Eq. ^ into the previous equation 
we obtain 



(4) 



where = 1- The expression of Eq. (|4]) is the 

most general formula we obtain for the heating rate and 
must be evaluated for each separate structure. In par- 
ticular, the spectrum of the fluctuating force has been 
obtained in Appendix \K\ as 



TreoC 



rKh + 1) 



v(/x-Im V X G(rA,rA,w) X V • /x) V , (5) 

where nth — Xj [e!^^ I^^t _ ^ jg |.]^g mean thermal pho- 
ton number and G(r^,r^,w) is the Green function. The 
Green function accounts for the spectrum of the fluctuat- 
ing magnetic field according to the following expression 



5B(r,r',c^) = (B(r,a;)B(r',c^')> 
= ^ n ("-th + l)Im 



(6) 



V X G(r,r',a;) X V 8{uj - J) 



All the information regarding the trapping parameters 
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and the geometry of the system are contained in the 
Green function. 



III. Heating rate 

We present our numerical evaluations for the heating 
rate of Eq. ([4]) for two different chip structures as shown 
in Fig.[TI^a)-(b). Both structures contain a superconduct- 
ing layer, either niobium or YBCO, in the Meissner state. 
We evaluate the heating rate for a ^^Rb atom held at a 
distance z above the chip. We consider the following 
trapping frequencies along the three directions in space: 
^v,x = 0.1 KHz and io^^y — lOv.z ~ 1 KHz. 

The heating rate along the jth trapping direction for 
the niobium-based chip structure can be obtained from 
Eq. ([4]) by substituting the Green functions for isotropic 
media as reported in [2a |. leading to 




100nm ^ iivn 10 ^/m 

Surface distance z 

FIG. 3: Heating rate for a Nb-superconducting chip at 4.2 
K, the thickness of the Nb layer is 1 ^m. Each of the three 
curves represents a different structure: the bare niobium chip 
(solid line), the niobium-|-gold chip with gold thickness of 50 
nm (dashed line) and a simple gold surface (dotted line). 



a, 
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where ax = cty = 1 and = 2, and /is is the Bohr 
magneton. The Fresnel coefficients R^^ and R^'^ , for 
TE and TM waves, are computed according to the na- 
ture and number of layers constituting the structure. In 
particular, the Fresnel coefficients for the nth layer follow 
the recursive formula 



R 



'm,m+l 



(8) 



where m = n -I- 1 denotes the consecutive layer to n, 
with thickness dm, Pm = i/fc^ — , = CmW^/c^, and 
Cm is the relative permittivity of the mth layer. The 
label q denotes either the TE or the TM components of 
the electromagnetic field. The ^ terms represent the 
Fresnel coefficients for the simplest geometry, two half 
spaces, and are given by 



I3n - (3n 



TM 



^mf^n "f ^nPn 



(9) 



In Figure [3] we plot the heating rate F,:^^ along the z 
direction versus the surface distance for three different 
chip structures at liquid helium temperature T = 4.2 K. 
We scale the heating rate by the Boltzmann constant in 
order to facilitate comparison with relevant energy scales, 
however this should not be confused with the thermody- 
namic temperature. 

The heating rate is highest above the gold substrate, 
while the bare niobium surface yields the lowest rate. 
At distances of 1 /im, we estimate the heating rate for 
the gold substrate to be of the order of F^^au — 10~^ 
K/s which is a measurable effect. At the same distance. 



the heating for the niobium-based chip is reduced but in 
principle still detectable with a top gold layer of 50 nm 
as Fe^Nb+Au — 10~^. However, the rate for the bare nio- 
bium chip is too small to be of experimental relevance 
because Fe^Nb — 10"^^ K/s. Trapping distances below 
1 /zm are challenging because of the van der Waals at- 
tractive potential. However, for z ~ 100 nm, the heating 
rate of the niobium-|-gold surface approaches the heating 
rate of the gold substrate. This indicates that at very 
small distances only the effect of the metal is relevant. 

For the YBCO-based chip structure, the Green func- 
tion for anisotropic media as reported in [2^ is adopted 
and the heating rate along the z-direction (perpendicular 
to the chip surface) is given by 



Meoc' 

„2 



(nth + 1) 



(10) 



with kl ~ Lu^/c^, while the heating rate Fej| for the plane 
parallel to the chip surface is given by the sum of the 
following expressions 



F, 



(nth + 1) 



(11) 
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drj 

in 



-2"^Im [B},\kl + v'Bl},] , 



with k1 . = 



— uj^^/c^. The scattering coefficients Blj-^^ 
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and pl^^Y the same role of the Fresnel coefficients of 
Eq. ^ but account for the anisotropic properties of the 
YBCO layer. These coefficients can be obtained by fol- 
lowing [23, ii]. 



X 



100nm ^ iim lOyum 100//m 
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FIG. 4: Heating rate for a YBCO-based chip at 77 K, the 
thickness of the YBCO layer is 1 /im. Each of the three 
curves represents a different structure: the bare YBCO chip 
(solid line), the YBCO+gold chip with gold thickness of 50 
nm (dashed line) and a simple gold surface (dotted line). 



Our simulations for the heating rate of the YBCO chip 
are presented in Fig. [H We plot the heating rate along 
only the z trapping direction because this term domi- 
nates. We observe that the bare YBCO surface yields the 
lowest heating rate, compared to the other two structures 
with the gold layer. However, the difference between the 
bare YBCO surface and the YBCO+gold chip structure 
is dramatically reduced at liquid nitrogen temperature 
T = 77 K, and is less than one order of magnitude for 
a wide range of surface distances. The heating rate is 
detectable for distances around 1 fim where Fe^ybco — 
re,YBCO-i-Au — 10~^ K/s. At Smaller temperatures, the 
rate corresponding to the two structures differs by a few 
orders of magnitude, i. e. Fe_YBCO+Au/Fe,YBCO — 10"^ 
at T = 4.2 K given by the fact that the YBCO pene- 
tration depth and the skin depth of gold scale with tem- 
perature. This results in the near field spectrum being 
dominated by the gold. As a consequence, at very low 
temperatures the presence of the gold leads to a marked 
increase of near field noise. This increase is less pro- 
nounced at temperatures closer to the YBCO transition 
temperature. We conclude that the heating rate close 
to a bare YBCO surface and to a YBCO-f-gold chip for 
typical trapping distances of a few /im does not result in 
a measurable difference at 77 K. Moreover, the typical 
values of magnetic bias fields and temperatures used for 
superconducting chip experiments suggest that a high- 
temperature superconductor is likely to be operated in 
the Shubnikov phase. This involves the partial penetra- 
tion of magnetic field in the form of vortices but this is 
not considered in the present work. 



IV. Spin flip lifetime 

In the following section the comparison between the 
bare superconducting chip and the superconductor -|-gold 
structure is extended by considering the spin flip lifetime 
r of a neutral atom. The transition rate for thermally 
induced spin flips Tsf = I/t between an initial hyperfine 
state \i) and a final hyperfine state |/) is given in its most 
general formulation as [2l[ 

Tsf = ^^^{f\S,\t){i\SM)SB{rA,u:sF) , (13) 

where ujgp is the spin transition frequency and {i\Sp\f) 
the spin matrix element for \i) — > |/). We restrict the 
calculation to a two-level system evaluating the spin flip 
lifetime for the transition \F — l,mF ~ —1) \F = 
l,mF = 0). We choose a typical experimental transition 
frequency losf = 27r560 KHz and we note that the 
near field noise spectrum does not change significantly 
in the rf frequency range. We evaluate the spin matrix 
elements via the Clebsch-Gordan coefficients and obtain 
for the non-zero matrix elements Kjl-S'a;!/)! — \{i\Sz\f) \ — 
1/8 (assuming that y is the spin quantization direction). 

The spin flip lifetime close to the niobium-based chip is 
obtained as the inverse of the following spin flip rate [2g| 

^^^•^'^ = 32A^ J {2nY 2 ]("th + l), 

(14) 

where R^^ is the Fresnel coefficient as in Eq. The 
spin flip lifetime as a function of the surface distance 
is plotted in Fig. [5] for three different structures. The 
shortest spin flip lifetime is obtained for the gold surface. 
The longest lifetime is found close to the bare niobium 
surface and is two orders of magnitude longer than the 
lifetime near the niobium+gold chip. 




lOOnm 1 //m 10 //m 100/im 



Surface distance z 

FIG. 5: Spin flip lifetime for a niobium-based chip at 4.2 K, 
the thickness of the Nb layer is 1 /im. Each of the three curves 
represents the lifetime corresponding to a different structure: 
the bare niobium chip (solid line), the niobium-|-gold chip 
with gold thickness of 50 nm (dashed line) and a simple gold 
surface (dotted line). 
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The spin flip rate above the YBCO-based chip can be 
written as j28| 



SF, YBCO 



{nth + 1) 



(15) 



drj ■ 



-2riz 



-Im 



B 



where Blj^ and Bjj2 the scattering coefiicients of 
Eqs. (|10lll2p . The spin flip hfetime t = l/rsF,YBCO is 
plotted in Fig. [6] for the three different structures: gold 
surface, YBCO surface, YBCO+gold structure. Simi- 
larly to the niobium-based chip structures, the spin flip 
lifetime close to the YBCO-|-gold structure is shorter 
than the one obtained near the bare superconductor, 
however the difference is less than an order of magni- 
tude. A complete accounting of the total lifetime must 
include other experimental restrictions, such as collisions 
with background gas Q . 
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100nm \ iim 10 //m 100/im 

Surface distance z 

FIG. 6: Spin flip lifetime for a YBCO-based chip at 77 K, the 
thickness of the YBCO layer is 1 /^m. The two upper lines 
represent the case with a bare YBCO wire (solid line) and 
with a gold layer of thickness 50 nm on top of YBCO (dashed 
line). The lifetime for a simple gold surface is reported for 
comparison (dotted hne). 



increasing gold thickness until it approaches the lifetime 
obtained close to the single gold layer. This happens for 
a layer thickness of the order of 10 /im. 
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FIG. 7: Heating rate and spin flip lifetime of an atom 1 [xm 
above a superconducting chip with a top gold layer. Both 
quantities are plotted as a function of the gold thickness d\ 
(as in Fig. [T] (b)), for a niobium-|-gold chip (solid curve), 
YBCO-f-gold chip (dashed curve) and for a gold substrate 
(dotted line). Both chips have a temperature of 4.2 K and 
the thickness of the superconducting layer is 1 pm. 



Conclusions 



We conclude this section by plotting in Fig.[7]the heat- 
ing rate and spin flip lifetime of the niobium-|-gold and 
YBCO-|-gold chip as a function of the gold thickness. The 
niobium-|-gold chip has the lowest heating rate. The rate 
increases for both structures with increasing the thick- 
ness of the gold until the layer thickness reaches 1 [ivn. 
The presence of the superconductor is then irrelevant 
and the heating rate is the same as what one would ob- 
tain close to a gold substrate. Typical values for the 
gold thickness are around 50 — 100 nm. For such values 
the heating rate for the YBCO-|-gold chip is one order 
of magnitude larger than the niobium-|-gold chip at 4.2 
K. Similarly, the niobium-t-gold chip exhibits the longest 
lifetime. For both structures the lifetime decreases with 



We have investigated the heating rate and spin flip 
lifetime for a superconducting atom chip considering the 
cases with and without a gold layer deposited above the 
superconductor. All the results presented are valid for a 
superconductor in the Meissner state. The heating rate 
for trapping distances around 1 /im is measurable and 
may become significant for a sufficiently cold atom cloud. 
However, by adding a top gold layer, this rate increases 
such that for layers thicker than 1 /xm, the presence of 
the superconductor becomes irrelevant. Even for a sub- 
micron gold layer, an atom held at a close distance of 100 
nm is influenced only by the metal, and both the heating 
rate and spin flip lifetime approach the values obtained 
in the proximity of a inflnitely thick gold substrate. 
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In particular, we have considered a niobium-based and 
a YBCO-based chip. At a chip temperature of 4.2 K, the 
difference between the case with the gold layer and with- 
out is marked both for heating rate and spin flip lifetime. 
Such difference is less pronounced when the YBCO chip 
is operated at 77 K. In summary, our study shows that 
the deposition of a metal layer above a superconducting 
layer can diminish the advantages of using a supercon- 
ductor to overcome the limitations given by a metal. The 
presented results are of interest for future development of 
technologies involving superconducting atom chips. Fur- 
ther investigations will include the treatment of type II 
superconductors in the Shubnikov phase. 



is the positive- frequency part, and the single component 
is given by 



€kjsdjGsn{r, r', w)/„(r', w). (A2) 



In the Heisenberg picture, the equation of motion for 
creation operator Oj' is given by 



k]{t) = - -^d.fikB, (rA) , (A3) 

and the bosonic field operator becomes 
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A. Quantum mechanical derivation of the heating 

rate 

We outline the quantum mechanical derivation of the 
heating rate based on the model introduced in Sec. [TTl 

The total Hamiltonian of a harmonically trapped atom 
held close to a substrate (the conducting properties of the 
substrate are not yet specified at this stage), is given by 
the sum of three different Hamiltonians as ~ Hf + 
Ha + Hint ■ The Hamiltonian for the atom in the unper- 
turbed harmonic trap reads Ha — fu^v,j + 5^ • The 
Hamiltonian of the electromagnetic field arising from an 
absorbing and dispersing medium is given within the for- 
malism of quantum electrodynamics for dielectric media 
[Hil as iJ/ = / cPrdujhujf'^{Y,uj)f{Y,uj), with f(r,w) 
and f'^{v,u}) being the bosonic operators accounting for 
the collective excitation of the medium and the electro- 
magnetic field. They satisfy the usual equal-time com- 
mutation relations f (r, a;), ft(r', (jj') —5{t — y')5{ijJ — lo') 
and their correlation function at temperature T reads 
(f(r,c^)ft(r',w')> = {nt^ + l)6{Y-Y')5{uj-uj'), where nth 
is the mean thermal photon number. 

The interaction Hamiltonian of Eq. ([T]) can be written 
in the rotating wave approximation as 



H^nt = - y (a]djiikB+{Tj^) + H.c, 



(Al) 



where Xj denotes either Xj^+ for n 



n+1 or Xj - for n —> 

1. The magnetic field i3^(rA) can be obtained via 

00 

B(r) = B+(r)+B-(r) where B+(-)(r) = / dwB(t)(r,w) 





/„(r, w,i) = -iufn 



(A4) 



-J "^^^ ^V<^ii^^^)dj £kqsdq G*„ (r^ , r , w) , 

where the time integral in the Markov approximation be- 
comes 



/n(l", t) = /„Jree(l", t) + i%(i)C(Wt,j - Uj) 



(A5) 



Xj fik UJ 

2 y/hneo 



^yei{r,uj)d.jekqsdqG*JrA,r, u) , 



and C(x) — Tr5{x)+iVx^^ {V denotes the principal part). 
This formal solution is going to be substituted into the 
magnetic field of Eq. (|A2p and making use of the follow- 
ing integral relation for the Green function [2^1 

/ d^s— e7(s,tj)G(r,s,t^)G*(r',s,cj) =ImG(r,r',w), 

(A6) 

we obtain 



B^ {va , W) = Bi. free (l"A , t^) + ia.j (t) 



1 



2 TreoC^ 



(A7) 



C(aJi,j — a;)9jlm V x G(r^,r^,w) x V 



ih 



[Ik 



where the gradient V = X]j=a; y z '^3 acting on the 
second argument of the Green function as denoted by 
the left pointing arrow. After integrating Eq. (|A7|) over 
bj and substituting it into Eq. (|A3p . we can write the 
equation of motion for the creation operator as 



a](t) = [-7j/2 + i{uJv,j - 5uj^j)] a] 
~Y2dj[ipB-j^^^{vA). 



(AS) 



where 7^ = 7j,n^n±i is the rate associated with the n — > 
nil transition along the jth direction. The rate arises 
from the 5{x) of the C,{x) function appearing in Eq. (|A5|) . 
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and is given by 

7j 



V (/J, ■ Ini 



(A9) 



At 



while the term Sujj of Eq. (|A8|) arises from the principal- 
part of the C(a;) function and denotes the frequency shift 
which we assume to be negligible. 

If the dielectric body is in thermal equilibrium with 
its surroundings, then the magnetic field is in a thermal 
state with temperature T, and the magnetic field spec- 
trum has to be multiplied by the mean thermal photon 
number nth + 1 • The spectrum of the magnetic field can 
be expressed as 



and hence the spectrum of the fluctuating force causing 
the shift of the center of mass is 



;inth + 1) 



(All) 



V(^/i-Im V X G(rA,rA, w) X V • /x^ V , 

such that the transition rate for the jth trapping direc- 
tion can be written as 



nx 



^j,n — *n±l 



2K^ 



SF{rA,i^v,j)- 



(A12) 



SB{rA,i^) 



;{nth + 1) 



Treoc^ 

Im V X G(r^, r^, w) X V , 



(AlO) 
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